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Abstract 

The response function of a homogeneous and dense hadronic system to a time- 
dependent (baryon) vector potential is discussed for holographic dense QCD (D4/D8 
embedding) both in the confined and deconfined phases. Confined holographic QCD 
is an uncompressible and static baryonic insulator at large Nc and large A, with a 
gapped vector spectrum and a massless pion. Deconfined holographic QCD is a diffusive 
conductor with restored chiral symmetry and a gapped transverse baryonic current. 
Similarly, dense D3/D7 is diffusive for any non-zero temperature at large Nc and large 
A. At zero temperature dense D3/D7 exhibits a baryonic longitudinal visco-elastic 
mode with a first sound speed 1 / \/3 and a small width due to a shear viscosity to baryon 
ratio rj/nB = h/A. This mode is turned diffusive by arbitrarily small temperatures, a 
hallmark of holography. 
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1 Introduction 



Solids respond to external stress elastically through their bulk and shear modulii K and 
// respectively, with almost zero dissipation. Liquids on the other hand, follow the lore of 
hydrodynamics with bulk and shear viscosities ^ and 77 accounting for dissipation. In contrast 
to the solid, the shear modulus vanishes in the liquid. The bulk modulus does not. 

This remarkable difference between solid and liquid disappears when the stress is time- 
dependent. Indeed, for a stress of finite frequency uj a liquid has a non-zero shear modulus 
much like the solid. In the long- wavelength limit, the dual description of a solid or a liquid 
follows from the visco-elastic equations with complex and frequency dependent elastic con- 
stants as we detail below. In this paper we will explore some of these ideas in the context of 
the AdS / CFT correspondence by analyzing the baryonic response functions at finite density 
for both D4/D8 and D3/D7 embeddings. 

The AdS/CFT approach [1] provides a useful framework for discussing large A^^^ gauge 
theories at strong coupling A = g'^Nf.. The model suggested by Sakai and Sugimoto (SS) 
[2] offers a specific holographic realization of hQCD that includes Nf flavors and is chiral. 
For Nf -C A^c, chiral QCD is obtained as a gravity dual to Nf D8-D8 brancs embedded 
into a D4: background in 10 dimensions where supersymmetry is broken by the Kaluza- 
Klein (KK) mechanism. The SS model yields a holographic description of hadrons in the 
vacuum [2, 3, 4, 5, 6, 7], at finite temperature [8] and finite baryon density [9, 10, 11, 12]. 

Hot and dense hadronic matter in QCD is difficult to track from first principles in current 
lattice simulations owing to the sign problem. In large Nc QCD baryons are solitons and a 
dense matter description using Skyrme's chiral model [13] was originally suggested by Skyrme 
and others [14]. At large N^ and low density matter consisting of solitons crystallizes as the 
ratio of potential to kinetic energy T = V / K N^ /p'^p ^ lis much larger than 1. The crystal 
melts at sufficiently high density with F ^ N^/Pp ~ 1, or sufficiently high temperature with 
F fa Nc/T fa 1. QCD matter at large Nc was recently revisited in [15]. 

The many-soliton problem can be simplified in the crystal limit by first considering all 
solitons to be the same and second by reducing the crystal to a single cell with boundary 
conditions much like the Wigner-Seitz approximation in the theory of solids. A natural way 
to describe the crystal topology is through with periodic boundary conditions. A much 
simpler and analytically tractable approximation consists of treating each Wigner-Seitz cell 
as with no boundary condition involved. The result is dense Skyrmion matter on [16]. 

At low baryonic densities holographic QCD is a crystal of instantons with the Wigner- 
Seitz cell approximated hy S^. The pertinent instanton is defined on S^xR [11]. At moderate 
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densities chiral symmetry is restored on the average with an equation of state [11]. This 
homogenous (on the average) hquid-hke phase is strongly coupled and not emmenable to 
standard Fermi liquid analysis. 

In this paper, we would like to follow up on the transport properties in the homogeneous 
phase originally discussed in [10] using D4/D8 to contrast them with some recent studies 
in [17] using D3/D7. In section 2, we recall the bulk characteristics of the homogeneous phase 
in D4/D8 and suggest that it may be identified with a strongly coupled holographic liquid 
prior to the restoration of chiral symmetry. In section 3, we derive the general formulae 
for the holographic currents induced by an external baryonic field in the linear response 
approximation for both D4/D8 and D3/D7. In section 4 we show that the transverse baryonic 
current for cold D4/D8 is saturated by medium modified vector mesons in leading Nc in 
agreement with [10]. The bulk static conductivity is zero. Large Nc D4/D8 is an insulator. In 
section 5, we develop the quasi-normal mode approach for hot and dense D4/D8 and D3/D7, 
both of which are conductors at large N^. For completeness we also discuss cold D3/D7 
in light of a recent result [17]. In section 6, we suggest a unified visco-elastic framework 
for interpreting gapless excitations in dense media in both the elastic (coUisionless) and 
hydrodynamic (collision) regimes. We argue that cold D3/D7 exhibits such a mode at large 
Nc with zero bulk viscosity and finite shear viscosity. In section 7, we suggest that the leading 
1/Nc correction to the baryonic currents in cold D4/D8 can be extracted from an effective 
baryonic theory using the Random Phase Approximation. Our conclusions and prospects 
are in section 8. A number of points pertaining to transport in dense holographic media are 
discussed in the Appendices. 



2 Homogeneous dense matter 

First we briefly review the bulk property of the homogeneous phase in the SS model. More 
details can be found in [10]. We consider Sakai-Sugimoto's original embedding [2, 3], where 
the D8-branes configuration in the r coordinate is constant and not affected by the existence 
of a background U{l)v field Aq. This corresponds to r = 6r/A, the maximal asymptotic 
separation between D8 and D8 branes. The DBI action of D8 branes with Aq is written as^ 

5dbi = -a I d^x j dZK^/^ ^ 1 _ bK^^dzKy , (1) 
^The integral is extended to (— oo, oo) to take into account D8 branes as well as D8 branes. 
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where 



a = 



6^2 



3%5 



b = 



3V 



(2) 



Mkk is the Kaluza- Klein mass and A is t'Hooft couphng. Now we introduce the baryon 
source coupled to Aq through the Chern-Simons term. We assume that baryons are uniformly 
distributed over space whose volume is V. For large A, the instanton size is 1/VX [4, 5]. It 
can be treated as a static delta function source at large Nc- For a uniform baryon distribution, 
the source is 

^source = N^Ub J d^X j dZ 8{Z)ko{Z) . (3) 

By varing the total action S'dbi + 'S'source we get the classical solution Aq: 



Ao(Z;ng)= / dZ 



nj2 



(4) 



which defines the baryon quark chemical potential fj,g as 



Hq{ng) = lim Ao{Z;nq) 

\Z\—>-oo 



(5) 



This relation also defines /ig as a function of Ug and vice versa. The baryon chemical potential 
fiB is 



ljLB^rnB + Ncjig , 



(6) 



where is the baryon rest mass. 

The interaction energy density eint, pressure P, grand potential Jl, and the baryon chem- 
ical potential ^b have been computed in [10], 



V 



dZK"^'-' 



a / dZK''/^ 



1 + _ 1 



4a26 



1 - 



jiB = ruB + Nc dZ 



f 

J — c 



A^c^b/4 



(7) 
(8) 
(9) 



^{abfR^ + bK^/^iNcUB/lf ' 
where V is the volume and Q is understood as a function oi hb through (5) and (6). At low 
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Figure 1: Sketch of the phases of cold D4/D8. 



densities they translate to 



2777^ N, 



liB ~ mB + 2lTx^^-^nB . (12) 
iV/A A^KK 

The baryonic contributions appear through the combination Nc/\Nf. The large Nc and 
large A limit are not compatible in the homogeneous phase. Compatibility with solitonic 
physics suggests that the large limit be taken first followed by the large A limit, which 
is also consistent with holography. This will be assumed throughout, unless specified other- 
wise. The homogeneous phase described by (8) breaks spontaneously chiral symmetry with 
density dependent vacuum-like modes [10]. In Fig. 1 we sketch the various phases of dense 
holographic matter at zero temperature. The low density part is inhomogeneous (solid) with 
spontaneously broken chiral symmetry, while the high density phase is homogeneous (gas) 
with restored chiral symmetry. Intermediate between the two is a possible liquid phase. 
Here we suggest that (8) may capture some aspects of the liquid phase still in the spon- 
taneously broken phase using holography. The solid phase binds with an energy density 
^int ~ dMnB where qm is the Madelung constant for the pertinent crystallization provided 
that the baryons are semiclassically quantized to account for the pion-interaction through 
the mesonic cloud [6]. The gas phase is homogeneous with eint ~ n^B^ and restored chiral 
symmetry [11] 
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2 . 1 Compressibility 

Holographic QCD at large Nc and large A is umcompressible. Indeed, under small scalar or 
longitudinal vector stress the baryonic density ns is expected to change locally to Ub + Sub 
so that the constitutive equations read 

MnB^=-Vp, (13) 
dtSuB + UBV -v^O , (14) 

by the Newtonian equation of motion (13) and baryon current conservation (14). The 
baryonic charges move with an acceleration [dP/dnB)/mB ~ 1/A which is suppressed at 
large A since trb = SnXNc [5]. Another way to say this is to note that (14) implies 
{df—cl'V^)6nB = 0, with the speed of the first or thermodynamic sound ci = •>/ dP/duB/mB- 
For the confined D4/D8 configuration 



21. UB f^^l 



1/2 



8 X^Nf J K I ^2 — ^ 



(15) 



after using (8). The bulk modulus is K = nsdP/dnB ~ n\{Nc/X)., with the compressibility 
X — 1/K {\/Nc)/n\. Holographic QCD is uncompressible at large N^. 



3 Holographic Baryonic Currents 

Baryon transport in confined D4/D8 occurs explicitly through 1/Nc effects. Contributions 
to the baryonic current to order A^^ are shown in Fig. 2. They follow from direct (a) or 
vector meson (b) such as the uj meson. All density effects in holography are suppressed 
at large Nc and large A. To illustrate these points, we streamline the dense analysis given 
in [10] using general notations to extend the results to finite temperature and also other 
brane embeddings. The induced metric on the D8 branes for both low (KK) and high 
temperatures(BH) can be written as 

c^^Ds = gttdt'^ + 9xxSijdx''dx^ + guudU'^ + gssd^l (16) 
= a (^^y\t' + (^^y\,,d^^^ (17) 
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q BjD 

^^^^)rW~^ OKJt^ ^^r^ 



q 



B h 



(a) 



(b) 



(c) 



Figure 2: Typical contributions to the baryonic response in D4/D8. (a) Direct A*"^ (b) Vector 
mesons A''^ and (c) Fermi (baryon) contributions A*""^ 



where for the KK background 



a 



KK 



u 



and for the BH background 



a^-/(f/), 7 



fiU) \dU 



fiU) 



Ut 

IT 



(19) 



The embedding information is only encoded in 7 and thereby guu- 

With the induced metric (16) and the pertinent guage fields, the general DBI action 
follows as 



^DBi = tr j d^xdU e-'^gls 



-googlxduu - glx^ouFou - glxguu ^ -Foi-Poi 



2 

'googxx 



FiuFiu — googxxguu ^ FijFij — Qxx ^ FijFijFuQFuo + 



i>j 



i>j 



1/2 

(20) 



where e~'^ = gs{U/ R)'^^^ and IN" = T8n4. The brane tension is Tg and is the volume 
of a unit S*^. ^ The F^ and terms cancel by symmetry. Among the F^ terms we only 
retained the relevant term for our discussion below. If we consider the fluctuation {Aa{x°')) 
around the classical configuration Aq (4), which is due to homogeneous matter at Z = 0, the 



^We absorb 27ra' into the gauge field for notational convenience. It will be recalled in the final physical 
quantities. 
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action can be expanded as 



f 



,(21) 



up to quadratic terms. F^fj = d^Ap — dpA^ — i[Aa) ^0] and 

1 



A = 



1 + 



(27ra')2 
07 



(A[,)^ 



(22) 



It is useful to change variable 



C/ = C/o(l + Z2)V3 ^ 



(23) 



where Uq is the coordinate of the tip of the D8-D8 cigar-shaped configuration in the confining 
background and the position of the horizon Ut in the black hole background. The range 
of Z is (0,00) contrary to U whose range is (C/0,00). Also this range can be extended to 
(—00, 00) if we consider D8 branes (—00, 0) together with D8 branes (0, 00) in a natural way. 
For convenience, we note the following useful relations 

+ U = U,K^I\ dC/=H|5_£_rfZ. /=l-(|)'i,. (24) 

where C/* is Ut for the black hole background. For the confining background, from here on 
and for simplicity, we follow Sakai and Sugimoto [2] and choose C/* = C^kk- In terms of Z 
the action reads 



S^NiY / d'^xdZki 



2AA'oF^o + A^FzoFzo + A/c2 FoiF^i 

i 



i>j 



(25) 



where we dropped the fluctuation independent part and the parameters {ki, k^jka, A, N) are 
different for each of the brane embeddings. They are summarized in Table 1. The case 
D3/D7 is separatly discussed below (section 5.3). We note that the dimensionless densities 



are d — 



ub for D4/D8con, d = WiS^^aW.'^B for D4/D8dec, and d = Jxn]n, '^q 



8 



D4/D8eon 

D4/D8dec 
D3/D7 



K 



N 

- 2167r3 



547r 

\NfNc 
2(27r)4 



h 

K 



k2 

X-4/3(27rr)-2 



h 
-1 

K 
-f 



d 



d 



d 



Table 1: Parameters of the different embeddings in (25). See text. 



zt 



D3/D7 with / = 1 — This exphcitly shows that the density effects are subleading at 
large A and large for fixed Nf. 

Now consider an abelian fluctuation in the Az = gauge 



A^ = a^{x\x\Z)+V^{x'',x') , 







(26) 



where vanishes at the boundary i.e. a^{x^, x^, oo) = 0, so that the boundary field is simply 
V^{x'^,x'^). Vi^{x'^,x'^) exists for all Z as the background. ^ With the Fourier decomposition 



(^fii^j ^^) ■ 
V,{Z,x\x') 



dwdq 
dujdq 



(27) 
(28) 



the quadratic action can be rewritten as ^ 



S^N l^dZ 



-\-(1t'^t(J'T — 2 Jt^tO'T — Jt^t^t 



where we introduced the gauge invariant variables 

Longitudinal mode : ol = qo-o + oja^, 
Transverse mode : ar = cuai, 



Vl = qVo + ujVs , 



(29) 



This is equivalent to the usual set up with the boundary condition A^{x , .t , oo) — Vp(a; , x ). The 
difference is in the equations of motion. The equation for is homogeneous but the equation for is 
inhomogeneous and sourced by as shown in (34). 

^For simplicity wc omitted the argument of the functions. Each quadratic term is a function of {uJ,q) 
(first) and {—to, —q) (second). We dropped the surface terms since they vanish on shell when the source is 
explicitly present. 
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with 02 = and used Gauss constraint A^a;aQ + k^qa'^ — 0. 02 = is a consistent choice 
since the transversal equation of motion decouples from the others. 
The differential operators Vl/t are defined as 

= ^^ A^^i^lf 2 ^^ + ^i^2A , (30) 
= -1 (dzkiksAdz - kM^^^ + ^~^hq^)) , (31) 
and the coefficient functions are 

_ kMAij^ + A-^ksf) 
JT = ^ ■ l^dj 

The equations of motion is 

DLttL = IlVl + 9L , I^TttT = ItVt. (34) 

With the formal solutions 

aL = 'Dl'ihVL + 9l) , ar = D?'(/tVt) , (35) 

the on-shell action reads 

S^-N J ^^dz(vLfL[Dl\fLVL)+VL]+VL[2fL'Dl'gL]+gL'Dl'gL 

+Vt/t[D?'(/tVt)+Vt]) . (36) 
The induced baryonic currents follow to leading order in large Nc and large A as 

Muj, q) = 2Nu; J dZfL[Dl\hVL + Ql) + ^l] , (37) 



Jt{uj, q) = 2Nuj J dZfT[DT\fTVT) + Vt] , (3^ 



where T)J^^ and T>^^ arc understood with the retarded prescription uj ^ uj + iO. fi, fr^QL 
are all recorded in (33). The longitudinal current involves Ql independently of Vj^ as Ql 
is triggered by the gradient of the baryonic profile Aq. This is the analogue of Pick's law 
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(baryonic charge diffusion). The terms involving Vl,t correspond to (Tl,t the longitudinal 
and transverse Fourier transforms of the space-time conductivities. The arguments (a;, g) 
are subsumed. 



4 D4/D8: Cold 

In the confined phase, the operators Dl,t are hermitian modulo the retarde prescription in 
frequency space. They can be diagonahzed using eigenmodes as discussed in [10]. Through- 
out the prescription a; ^ a; -|- iO is subsumed. 



4.1 Longitudinal Mode 

The longitudinal operator (Dl) is 



When 5 = or a; = it is easily diagonahzed, since 

D^(u; = 0) = -^dzKA^dz + K-^'^A . (40) 
T 

The Green's function (D^^) may be expanded in terms of the complete set of eigenvalues 
that diagonalize 

^L{q = 0)/ = [—^) A/ , 

^l{u; = 0)/ = [-^) A/ , (41) 
where ^ and ^ are weight factors. Using the complete sets, 

{dzKMz) Xn = -(i^-^/^A) {XlfXn , 

(dzKA'dz) = -(K-'/'A) , (42) 
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Vector modes 



Pion and Axial vector modes 
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Figure 3: Vector mode functions (Left) and axial-vector and pion mode functions (Right) 



we have 



{Z\T),\q = m') = Y.^ 



^2 + {Xlf 



neN 



which are the results of [10]. Typical behaviors of Xn and ^„ are shown in Fig. 3. 
For small ij,q we may write without loss of generality, 

(Pn{Z)(pn{Z') 



ne2N 



ne2N+l 



+ c2g2 + (AX)2 



(43) 
(44) 



(45) 



The first contribution is from the density dependent axial-vector mode, the second contri- 
bution is from the density dependent pion mode (strictly speaking its U(l) partner at large 
Nc), and the last contribution is from the density dependent vector mode. The denominators 
are the dispersive modes, while the numerators capture their residues. The even-odd in the 
labelling of the modes translates into odd-even in the parity of V9„(Z). The baryonic current 
reads 



Jl = 2Nu j dZK~'/-'A [D-^\K-'/-'AVl + 9l) + Vl] 



(46) 
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The Z-integration picks only the vector or even- modes of (43) since Vl is trivially even. The 
longitudinal baryonic conductivity in the confined case is 

aL = 'iNu J dZX-V^A [1 + D^^X-i/^A] . (47) 

The longitudinal mesonic propagator admits the mode decomposition (45). Prom (42) 
it follows that 

/ dZhxn = ^ fdZ (dzKAdz) Xn = ^ {KAdz XntZ = 

J n J n 

is a zero boundary term. The longitudinal baryonic conductivity simplifies 

aL = -2Nu; J dZf^ , (48) 

and so does the longitudinal current. The longitudinal conductivity vanishes at a; = 0. 
Confined holographic QCD is a static insulator at large and large A in agreement with 
our recent analysis [12]. 
We now note that 



= aI ' "'^ = :f " V JdZK-^A-^' ^^^^ 

where and /J have been derived in [10]. At high density the pion speed vanishes as 
~ I/tib- The propagation of the axial charge stalls in very dense matter. For small 
momenta q the poles develop at 



while for small frequencies a; 

4(-(Ai)^) + (A^r = =^^n= (f )' ' (51) 

since = — (A|)^ from (44). In the confined D4/D8 embedding, the vector and axial modes 
disperse through 

uJn^Xi + -^^q\ (52) 
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Figure 4: Dispersion relation for vectors (Left) and axial- vectors including the massless pion 
(Right), n = 1.26no 



where is the rest mass and ^^tt- is the kinetic mass (Fig.4). To this order, the imaginary 
parts vanish in holographic QCD [10]. Indeed, vector, axial-vector and pionic modes are 
expected to be absorbed by excited and/or recoiling baryons which are l/Nc suppressed 
effects in cold and dense QCD. 

4.2 Transverse Mode 

The transverse operator {Dt) is given by 



'Dt = ^ {dzKAdz + K-^'\Auo^ - A-ig2)) 



(53) 



For g = it diagonalizes trivially through 



1 



Vriq = 0) = :^dzK/\dz - K-^'^/\ 



(54) 



The Green's function (2)^^^) may be expanded in terms of the complete set of eigenvalues 
that diagonalize 



DT{q = 0)/ = - ( ) A/ 



(55) 



by using the eigenvalues 



{dzKMz)Xn = -{K-'/'A) {XlYxn . 



(56) 
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The Green's function is then expanded as 

= Q)|g) = J2 ^gl'^'lf ? + , (57) 

which is the same as Dj^^. A rerun of the arguments for the longitudinal current response 
yields 

Jt{u, q) = 2Nu J dZMV^'ifrVT) + Vt] , (58) 
where again the retarded prescription is subsumed and 

In the confined phase, the transverse conductivity follows 

(7t = 2Nuj J dZfril + 1»^Vt] ■ (60) 

Using the vector meson mode decomposition for D^^ (45) and the relation (151) we can 
simpUfy the transverse conductivity 

(7T = (7L + 2Nu; j dZ{fT + /l)[1 + ©^'(/r + h)] , (61) 

with 

fT + h^-^.K-y^A-K (62) 



a;2 



The transverse conductivity is vector meson mediated as shown in Fig. 2. For q — 0, ax — ctl 
and vanishes for a; — > in agreement with [12]. 

5 Quasi-normal mode analysis 

we now turn our attention to the deconfined phase of dense holographic models with non- 
hermitean or absorptive boundary conditions. For that, the retarded prescription on the 
inversion of is best captured by the quasi- normal mode analysis. The latter is enforced 
analytically by matching for the gapless modes and numerically for the gapped modes. We 
now present the general formulas pertinent to the longitudinal and transverse currents. 
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For u <^ l,q <^ 1 the equations (30) and (31) are reduced to 



{dzhksAdz) ot = . 



(63) 
(64) 



The general solutions are 



poo / 2 2 

= aL{Z){q 0,a; ^ 1) , 



(65) 
(66) 



where we imposed the vanishing Dirichlet boundary condition at the boundary Z — oo 
i.e. aL/T(oo) = 0. Cl/t will be determined by imposing incoming boundary condition at 
Z = which corresponds to the location of matter (confined phase) or the black hole horizon 

(deconfined phase). 

To constrain Cl/t we need to know the behavior of Ol/t around Z = 0. First, wc solve 
the equations (30) and (31) near Z = with fixed uj and q. Second, we take the limit 
a; <C 1, 5 -C 1. A^aL = and AtOt = may be written as 



ttr + 



A'(3h(f + ^'2A2(l - A-,(1/A-,)'A/A0) K 



A{ksq^ + u;^A'') 



, , LU^A- + k-iq^ , , 

- 1^2 JT^^^ ttL = ,(67) 



k,A^ 



{kik.AY 
kiksA 



LO^A^ + k^q^ 



Or = . 



(68) 



In this analysis the variable Z introduced in (23) is not convenient due to the complicated 
form of A. To make it simpler we introduce the new variable: 



1 



z — 



;i + z2)V3' 



(69) 



which is nothing but Uq/U in terms of the original coordinate in (17). In this coordinate the 
boundary is m = and the horizon or the matter location is m = 1. 



5.1 D4/D8: Cold and Dense 

Before analyzing the deconfined phase of D4/D8 with black-hole background absorption, it is 
amusing to ask whether the D4/D8 confined background with matter in the KK background 
could be also addressed with absorptive or non-hermitean boundary conditions. After all cold 
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matter disperses and absorbs waves much like a black-hole. Prom (65) and (66) it follows 
that 



(70) 
(71) 



where 



a^{z) = / dz' 



1 



1 + d^z'^ 



/3 



(72) 

(73) 



for u,q <^ 1. Recall that the density d is defined by ^2^^\fS below (25). At small 
densities we expand 



aL{z) = Cl 
with 



r/2 Q^4 

[iv^ - q')a^'\z) - -{u^ - 3g2)a(i)(^) + —{u^ - 5q')a^'\z) 
2 8 



+ 0{d^) , (74) 



























1 



















dz' 
dz' 



Z' 



1-z'^ 3 
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^arcsin(\/l — z^) — — ^ , 



2;' 



1-z'^ 



z 



120 r(2/3) 



z' 



JW 



1-z'^ 



We don't show a'^'^^z) cxpHcitly since it is long and not illuminating. 

To impose the incoming boundary condition we expand the solution around z = 1 



aL{z) = ClA + ClBVT^ + 0((1 - zf/^) 



(75) 
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where 



3 



B 



240r(2/3) 

-8 + 12(f - 15d%^ + (8 - A(f + 3d^)cu^ 



35841(4/3) 



(76) 



On the other hand we may first solve the equation near z—1. In general A^ql — and 
AxttT — are 



1 



2{1 + d^z^){q^ -u\l + d^z^)} 2(1-^3) 2z 



3^2 



(5d^z^) 

2(1 + ^2^5) 2(1 - z^) 2z 



9{g"-a;"(l + dV)} 

4z(l -Z3) (1+^2^5) 

, 9{g2-a;2(l + dV)} 

4^(1-^3) (1 + ^2^5) 



= . (77) 
= , (78) 



and reduce to 



1 1 ^30^^-^ 



near z = 1. The solutions are 



(79) 



aL/T = C,e-V37(i-^) ^ 



_ 2 T 
7 = 0; — 



(80) 



1 + ^2 ' 

where we imposed the incoming boundary condition aX z = 1. Their expanded form reads 



aL{z) ^Ci- %Ci^\fY^z + 0(1 - z) , 
^ = uj^-q^ + q^d^-q^d^ + 0{d^) . 



(81) 
(82) 



A comparison of (75) and (81) yields 



B + iAy/z^ = , 



(83) 



and the dispersion relation is 



a; = ± 



q + 0{d^)q' + 0{d'')q-^ + 0{q^) . 



(84) 
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The latter resums to 



(85) 



This is consistent with the zero mode result obtained in (45), (49) and Fig.4 where we also 
found the zero mode solution odd in Z. Interestingly enough, the quasinormal mode analysis 
when applied to the confined and dense KK background with absorptive boundary condition. 



it yields a massless pole which is the pion pole with a speed c-,^ — 1/ \/T+~d?. Note that there 
is no imaginary part. The reason can be traced back to the + (outgoing) and — (incoming) 
wave assignment in (83), both of which solve 



= S ± iA^/3-f 



- 4^2 + 3^4 
4vl 



^/J±iV3 



x/7 = 



^6) 



for 7 = 0. 

For the transverse mode we may follow the same procedure with the substitution in (75) 



ariz) = aL{z){q ^ ,a; ^ 1) . 



(87) 



The relation (83) yields 



tt ^ 407r3/2r(2/3) j ^ + > 



+ 0{q) , 



(88) 



which shows that there is no massless excitation. This channel is indeed gapped in the 
confined D4/D8 case as we discussed earlier for the case of refiecting boundary conditions. 



5.2 D4/D8: Hot and Dense 

The absorptive boundary condition is more appropriate for the deconfined BH background 
that we now discuss. For that, we rerun the same steps as we did in the previous section. 
First we solve the equations for tn <^ 1, q <^ 1, where to = and q = From (65) and 
(66) 

aL = Ci(ro'a„(^) - q'aq(^)) , (89) 
ar = CraM , (90) 
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where 



dz' \ 



10 

2 

15' 



,3/2 



(1 + d2^'5)3 

3 



(91) 



+ 2 2i^i[3/10, 1/2; 13/10; -dV^ 



73 



(92) 



Recall that the density d is defined by -^iFxi^W/TW^^ below (25). To impose the incoming 
boundary condition we expand aL around the horizon. 



0^(1 - e) = aL{l) - ea'^(l) H 

= CLXo'a^{l) - ^^=f - CLq'a,(l) + 0(6) , 

where 0^,(1) has a logarithmic divergence and aq(l) is finite. 
In general A^ol — and AtOt — are 



(5d22'^){3(l - ^3)^2 _ tp2(^ ^ ^2^5)1 
2(1 + d^Z^){{l - Z3)q2 - tt)2(l + (^2^5)} 



3-2 



(l-23)q2_tx,2(l + (^2^5)l_^3 2Z 

9{(1 - z3)q2 - tt)2(l + c^V)} 



4z(l-;z3)2(l + 0(2^5) 



-ol = , 



Ot + 



3^2 



1 



(93) 



(94) 



9{(l-r^)q2-m2(l + dV)} 

- 4.(1 -.3)2(1 + ^2,5) = ,(95) 



_2(l + d2.5) (1-.3) 2. 

which simplify to 

near the horizon with . = 1. The incoming wave solution near the horizon is 



(96) 



aL,T^{l-zr^/'F{z), 



(97) 
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with F{z) a regular function near z — lorz — 1 — e. Assuming tt) Ine <S 1 we have 



e-i^/^F(l-e) = F(l-e)-— lneF(l-e) + 

2 



lne|,^oi^(l) + 0(e) 



(98) 



By comparing the singular part of (93) with (98) we get 



(99) 



By comparing the regular part of (93) with (98) we get the dispersive retation for the 
longitudinal baryonic waves 



^ _ _iV0 _ i3aq(l)yrT^ q^ 
~ 2" 2 m' 

For small ro and q but fixed q^/ro the dispersion relation is 



(100) 



UJ 



2 2^15 LvTT^ 
g2 / 2^2 16^4 



27rr 



13 299 



+ 2 2Fi[3/10,l/2;13/10;-ci' 



(For small d) 



.j^ / 2r(i/5)r(i3/io) , r(i/5)r(i3/io) 

'2nT V 5r(l/2) ^ 5r(l/2) 



-8/5 



+ 



(For large d) , 



(101) 



where both the small and large baryon density limits are displayed explicitly. The longitu- 
dinal diffusion constant is 



Dj 



VT+d? 

27rT 



^^2|= + ^,Fi[3/10, 1/2; 13/10; -d^] 



(102) 



For zero baryon density this is D = 1/2'kT. In the dcconfined phase of D4/D8 the baryonic 
charge diffuses whatever the density. This is expected from baryon number conservation. 
The presence of the BH in the deconfined phase overwhelms the Fermi effects. 

A rerun of the analysis for the transverse baryonic current follows the substitution 



ar = aL(q ^ ,tt) ^ 1) 



(103) 
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Comparing the singular part of (93) and (98) gives 

c. = !5lHvTT5f(:), (104) 

and comparing the regular part of (93) and (98) yields 

tt)3 = i2 . (105) 

Thus there is no hydrodynamic pole. The transverse baryonic current in dense and 
deconfined D4/D8 is still gapped. It is much like a transverse plasmon. 

5.3 D3/D7: Hot and Dense 

For comparison, let us consider in this case the non-chiral and non-confining embedding with 
D3/D7 at finite temperature and finite density. We consider the massless quark embedding 
where analytic solutions are available [18]. The induced metric becomes simply AdSr, x 
independent of the gauge field. 

ds'^^^{-fde+dx')+r'^^dz^+R'dni, / = (106) 

where R = ATiggNca'"^ is the curvature radius. We work in units of -R = 1. Zh = ttT 
where T is the temperature. SUGRA and SYM quantities will be tied by a' = 1/ VX with 
A = AnQsNc- With this metric, we compute the DBI action as 



-S'dbi = — ^ tr d xdZ g^g 



'dmQxx ^ FizFiz — goodxxdzz ^ FijFij — g^x ^ FijFijFzoFzo + 

i i>j i>j 



1/2 

(107) 



The result is analogous to the D4/D8 case (20) with three differences: 1) 3\r = T^fl^; 2) there 
is no contribution from the dilaton; 3) g^^g appears instead of g^^g, since the compact space 
is not 51 

To consider finite baryon density (or chemical potential) we set the background vector 
C/(l) field Aq{Z) in bulk. Its form follows from minimizing the DBI action (107): 

2WA' = , ^ . (108) 



22 



(2nf 



-rig [18, 17]. 



We explicitly recalled 27ra' and d _ 

Following the analysis in D4/D8 above, we now consider mesonic fluctuations around the 
density background Aq. In the general form cast in (25) the quadratic action reads 



s = mT / d^xdzki 



i 

+Ak3 PizFiz + A-^k2h Y ^ii^^ 



i>j 



(109) 



where the information of the background field Aq is encoded in A and 



N 



XNfN^ 



A = Vl + d^Z-^ , 



2(27r)4 ' 

k,^z\ z-'r' , k' = f-\ 

as in Table. 1. 

In this general form we can use (63)-(68). In terms of the variable z 



and d = 



27rr 



we have 



ql = CL(m^a„(z) - q^aq(z)) , 



(110) 



27rr' 



(111) 

(112) 



with 



dz' 

VH- d 



Ottjl^;) = dz =- . 



(113) 



(114) 



Note that the integrand in a^, exhibits explicitly the BH horizon ai z' — 1 in units of 
temperature. At zero temperature the integrand smoothly reduces from 1/(1 — z'^) to 1. 
However, the BH singularity makes the integral logarithmically divergent at the horizon. As 
a result, the zero temperature limit is singular and will be considered separatly next. To 
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impose the incoming boundary condition we expand ql around the horizon. 

aL(l-e) = ai(l) - ea'^(l) + --- 

= CLtv'aM - ^^=^ - Cm'a.il) + 0(e) , 

Alternatively, the zero mode equation (Dl^l = 0) is 

Sd^^s /3(1 - z^)q2 - tv\l + d^z^) 



d.ttL + 



(1 + d2^6) \^ (1 _ ^4)q2 _ 1^2(1 + ^2^6) 



1" 



(l-2^)q2-ra2(l + d226)l_^4 ^_ 



+ 



4 



3d2^5 



1 - 1 - 1 + d2^6 

4^3 1 



ol = , 



(1 + d2;26) 1 _ ^4 ^ 

4 



+ 



1 - ^4 v 1 - z4 1 + d^^e 



which reduces to 



near the horizon at z = 1. The incoming wave solution is of the form 



jaL/T = , 



with F[z) a regular function near z = 1. For txi Ine -C 1 we have 

e-i^/2F{l-e) = F(l-e)-— lneF(l-e) + - 
= F(l)-^lneUoF(l) + 0(e) 

A comparison of the singular part of (115) and (120) yields 



tn 
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and a comparison of (115) and (120) yields 



1 = - i2a^{l)VlT^^ . (122) 



For small tt) and q with fixed q^/tt) the dispersion relation follows 



J \ 2 2 



^ ^ -M/l+lp^J ^2Fi[3/2,l/3;4/3;-d^ 

~ -i^r-T^ 1 + :^ , - TTT^ , H For small 



27rr I 8\{nTfJ 112 \{nTf J J \ (ttT) 

. f 2r(7/6)r(4/3) ( d \ ^'^ r(7/6)r(4/3) / d ^ '^'^ 

[ r(i/2) ivprJ + r(i/2) 

For Large ^ ) , (123) 



The longitudinal diffusion constant for hot and dense D3/D7 is 



1 L f d 



Dl - ^yi+^(^J 2F,[3/2,l/3;4/3;-d^ (124) 

As mentioned earlier the zero temperature limit is singular owing to the occurence of the 
BH pole in the issuing integrals. 

To analyze the transverse baryonic current ^ in the same limit of small o;, g we follow the 
same procedure with the substitution in (115) 

Or = aL{q ^ , tt) ^ 1) . (125) 
A comparison of the singular part of (115) and (120) yields 

Ct = i2roVl + d2F(l) , (126) 
while a comparison of the regular part of (115) and (120) gives 

tt)^ = i2 , (127) 
which is incompatible with the limits. The transverse baryonic mode in hot and dense D3/D7 



^See [19] for related work on the dispersion relation 
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is gapped much like the transverse plasmon in dense matter. This reflects on the long-range 
nature of the transverse forces in holography. We will comment further on this point below. 



5.4 D3/D7: Cold and Dense 

In a recent analysis [17] have reported the occurence of a zero soundm^ode in cold D3/D7. For 
completeness we now rederive their results using our general result (65). At zero temperature 
we set Zh = in (106) and change the variable to z — 1/Z. 
For a; <^ 1 and g <^ 1 (65) gives 



Ql = CLiuP'a^iz) - q^aq{z)) , 



(128) 
(129) 



with 



aq{z) = / dz' 
Jo 



-^==3 = l^A 2i^i[3/2, 1/3; 4/3; -z'd']) , 



lo VI + d^z'^ 2 
Near the horizon is expanded as 



dz'- ^ - .6^21 



= -^^(2Fi[l/2,l/3;4/3;-^V]) , 



(130) 



aL = CLA- + CLB + 0{l/z^) , 
z ^ ' 



with 



A = 



B = 



_ (g2_3^2)^-2/3r(i/3)r(i/6) 



d ' 18r(l/2) 
On the other hand the zero mode equations Dl/t(Il/t = are 



(131) 



(132) 



d^aL + 



■2^5 



1 



W^Z^ _ 1 

:i + ^2^6) - - 



dzC^T + { oj'^ — 



1 + d'z^ 



ar — ■ 



1 + d^z^ 



ol = , 



(133) 



Near the horizon the incoming solution is 



aL(z) = Ci- 



(134) 
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and for ojz 1 

dLiz) — \- icuCj , (135) 

z 

A comparison of (131) with (135) yields 

Aiuj^B , (136) 
which yields the dispersion relation reported in [17] 

for a massless excitation. Holographic D3/D7 at arbitrarily small temperatures is diffusive. 
It is not at strictly zero temperature with the occurence of a long-range collective mode. In 
the next section we suggest that this is a visco-elastic mode, and thereby generalize it to 
massive quarks. Any amount of temperature (collisions) destroy the Fermi-surface at large A 
and large Nc- Indeed, while the temperature effects are of order through the underlying 
BH metric, the density effects are 1/A and 1/Nc suppressed through the Np embeddings 
either D7 or D8. 

Finally and for completeness we note that the transverse baryonic mode follows also from 
(131) with the substitution 

aj. = aL{q ^ , u; ^ 1) . (138) 

From (136) we get 

r(i/2) ^ ^ 

The transverse baryonic current is gapped in cold D3/D7 much like the current is plasmon- 
gapped in a metal. 



6 Visco-Elastic Analysis 

The occurence of a collective mode in cold D3/D7 suggests that collectivity through the 
possible occurence of a Fermi surface at strong coupling maybe at work. To understand 
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Figure 5: Visco-elastic domain versus the free particle-hole continuum. See text. 



that, we propose to understand this collectivity by unifying the hydrodynamical or collision 
regime with the elastic or coUisionless regime. 

In Fig. 5 we show different propagating domains for a wave of frequency lo and momentum 
q in liquids. The dashed curves are typical wave dispersions. The free particle-hole continuum 
occupies the lower quadrant, r is a typical relaxation time to equilibrium, say r ~ l/T (hot) 
and T ^l/ji (cold) for conformal and strongly coupled theories. For waves with ujt ^ 1 and 
large velocities compared to the Fermi velocity vp of a quasiparticle, we expect coUisionless 
wave propagation or elastic regime. For waves with ujt but still large velocities compared 
to the Fermi velocity I'p we expect collision wave propagation or hydrodynamic regime. 
Typical cold media behave elastically at low temperature and hydrodynamically at higher 
temperature. Thus, an elastic mode can be turned inelastic by just raising the temperature. 
This is typically what happens in liquid He^ where the zero sound transmutes to the first or 
thermodynamic sound by changing its frequency or temperature to interpolate between the 
coUisionless and collision regimes. To understand these regimes we now introduce a unified 
visco-elastic framework. 



6.1 Generalities 

In a homogeneous and isotropic solid, the constitutive equations for the elastic displacement 
field u{t, x) are discussed in the canonical book by Landau and Lifshitz on the theory of 
elasticity [20]. Specifically 

niB = (X + (1 - 2/p) ii) V(V • u) + iiV^u + F{t, x) , (140) 
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where K and // are the bulk and shear moduh, p = 3 is the dimensionahty of space, F is an 
external volume force, is the baryon equilibrium density and is the bare baryon mass. 
The massless case more pertinent for D3/D7 will be deduced by inspection below. Without 
loss of generality, we may write 

F(t,f)=nB^^^, (141) 

— * 

which is the 'baryon electric field'. Since the transverse part of A induces a 'baryon magnetic 
field' we expect (141) to also include the magnetic contribution as it plays the role of the 
Lorentz force. Since wc are interested in the induced baryon current through A, the magnetic 
effects are second order and will be omitted. The baryon current density is 

J(t,x)=,.«^. (142) 

Inserting (141) and (142) in (140) and taking the Fourier transforms yield 

-iujmBj{uJ,q) = h i 



Li -; -> 

H — j{u},q) - iunBA{LJ,q) . (143) 

Decomposing the current j = jx + Jl and the potential A — Al + At along q and transverse 
to q allows the transverse and longitudinal induced currents 

M^, Q) = . ""^iZl^ M^, , (144) 



and 



Jl(^, q) = ^- 7^""^. , A,{uj, q) . (145) 



p I rig I niBijJ^ 



Jl relates directly to the induced baryon density through the local conservation law as qji^/uj. 

The response currents (144) and (145) have a direct analogy with their counterparts in a 
liquid. Indeed, using hydrodynamics for the baryon current density in a liquid we can write 
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the analogue of (143). In the hnear response approximation 

-iu;mBj{uJ,q) = h i-- : 

\nB ub \ pJ Ub J lOJ 

+ ^— — j{uj,q)-iojnBA{oj,q) , (146) 

Hb lUJ 

where the hydrostatic pressure term p in the Euler equation was traded with the longitudinal 
baryon current through the continuity equation, 

q) = -Q^VuBiu, q) = -q [-—^ J , (147) 

where 7] and C, arc the shear and bulk viscosities; p is the equilibrium pressure as a function 
of density with K = nBdp/dn the bulk modulus. (143) is very similar to (146) except for: 
1/ the shear modulus in the solid becomes purely imaginary or —iur] in the liquid; 2/ the 
bulk modulus acquires an imaginary part through —iuj( in the liquid. Both imaginary parts 
vanish at a; = making the liquid insensitive to shear at zero frequency. This also means 
that their contributions in j^,r are diffusive. 

The solid and liquid visco-clastic coefficients can be described in a unified manner through 
(144) and (145) by substituting 

K K(a;) = K{u;) - iu;C{uj) , (148) 
/X M{u;) = i^{u;) - iujr]{uj) , (149) 

as the complex and frequency dependent bulk K and shear M visco-elastic coefficients. The 
solid has real M with a small imaginary part rj (zero up to the uncertainty principle) while 
the liquid has imaginary M with a small /i and large imaginary part r). The bulk modulus 
K is about the same in liquid and solid, and of the same order of magnitude as the shear 
viscosity r). 

In hght of (149) it follows from (145) and (144) that the longitudinal current admits a 
gapless pole (compression mode) at 

. . (^ ^ 2 (l - 1) ^)"\ - i (^ ^ 2 (l - i) ^) f (150) 
yruBUB \ pJ mBUBj \mBnB \ pJ mBUB J 2 
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while the transverse current admits a gapless pole (shear mode) at 



q — i- 



rUBUB 



TJlBflB 2 



(151) 



We see that for finite frequency waves and in the long- wavelength approximation the way 
a solid responds to external wave-stress is similar to the way a liquid does. The difference 
is that in a solid /i K and (, r} are small, while in a liquid /x is close to zero and V are 
large. In the hquid the transverse or shear mode becomes diffusive. 

6.2 Cold D3/D7 

D3/D7 at finite density yields a gapped transverse baryonic current and a gapless longitudinal 
baryonic current [17]. The gapless longitudinal baryonic current can be compared with the 
longitudinal visco-elastic mode (150) with 



PJ rnsUB 



uibUb 

c 

TUbUb 



+ 2 



1-i 



1/2 



7] 



PJ TUbUb 



PI^B 



(152) 
(153) 



The compressibility is readily tied with the equation of state for any embedding 



K 



TTIbTIb 



dP 



1 
P 



(154) 



since e — pP = in a conformal theory. ^ The energy momentum tensor is still traceless at 
finite temperature and density for massless fermions. The gapless and longitudinal baryonic 
mode has the speed of the first sound ci = 1/ ^ for zero shear modulus /x = and massless 
quarks. For massive quarks e — pP ^ 0. For D3/D7 it follows from [18] that 



e = -7iV(2W)^(/x^ - m^)(3/i^ + m^) , 



(155) 
(156) 



is the dimensionality of space and P is pressure. 
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where 7 ^ 0.363. The visco-elastic mode has a speed 




c. = \ =W:fV-:^, (157) 



in agreement with the detailed quasi-normal mode analysis in [21]. 

Since the mode to = q/ ^/S lies within the free particle-hole continuum as shown in Fig.2, 
it is susceptible to Landau-like damping through single particle-hole or multi-particle-multi- 
hole. Again, the visco-elastic analysis suggests that the bulk viscosity in cold D3/D7 with 
massless quarks is 

^ = -r^ , (158) 

riB 2(p-l)' ^ ^ 

where we used (174) and (175). In conformal theories the shear modulus (/x = 0) and bulk 
viscosity are zero {(). For p = 3, we have tj/ub = h/i, where h has been restored. This is to 
be compared with h/ 6tt argued in [22] using the Stokes- Einstein relation and the uncertainty 
principle for cold and strongly coupled Fermions. 
It is worth noting that (158) can be rewritten as 

V _ f P \ h _ nphp h ^^^g^ 



VP ~ 1 / 2p NcNfUBP ' 
where 

AT AT f^P dPk 1 i^T 

np_iycJyfio j2^2k_N^Nf p 



ub fkp jjPk_ 2kF p-l ' 

Jo {2n)P ^ 



(160) 



is the ratio of the quark density at the Fermi surface normalized to the baryon density. For 
massless quarks UBkp = e + P with Uq = NcNjUb so that in general 

^ -^^^ (161) 



e-\- P Uq p 

This result can generalized to finite mass by noting that at zero temperature e + P = Ub/J'B 
and substituting 2k —>■ 2y^P~+m| in in (160). Specifically, 



2Fi(l/2,p/2;l+p/2;-4) , (162) 



ub ~ j^" ~ 4:k 
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with Vf — kF/rUq. Thus 

- = {^) 2^i(V2,p/2; 1 +P/2; -vl) ^ . (163) 
The visco-elastic (137) for massless quarks, turn to 

u;^c,q-'-(^l-^^ 2Fi(l/2,p/2;l+p/2;-4)^ ^ , (164) 

for massive quarks, ci is the first sound speed. For DS/D7 it is exphcitly given in (157), 
while for arbitrary Dp/ Dq it follows from the known equations of state [18]. 

In D3/D7 any infinitesimal tcmpcratm^c washes out the Fermi surface, resulting in a dif- 
fusive baryonic phase. Thermal collisions at strong coupling take over the collisions through 
the Fermi surface however small is the temperature. As noted earlier, this is the hallmark 
of strong coupling holography whereby the BH contribution is of order A'^^ while the Fermi 
contributions are 1/iVc suppressed. 



7 Random Phase Approximation 

If cold D3/D7 exhibits collectivity at large A and large that is consistent with a consti- 
tutive visco-elastic analysis, why the cold D4/D8 results above are all gapped. The short 
answer is that in D4/D8 the baryons are solitons, so baryonic motion with a Pauli- blocked 
Fermi surface is subleading in 1/Nc as shown in Fig. 2. Baryons move semi-classically by 
quantizing the isorotations and rotations both of which are 1/Nc suppressed. To estimate 
some of these contributions we note that the baryons in D4/D8 are flavored instantons with 
core sizes of order 1/-\/A. They are heavy with uib — Stt'^NcX. So the semiclassical descriptive 
of the translational zero modes follow from the point-like effective action 

S = <P^ (^td, - -I^B^<P- ^«(</.V)' + •••• , (165) 

where the repulsive interaction a = (247r^/4M|^^)(A^c/A) is set to reproduce the energy 
density and pressure of the holographic matter (8) with p ^ cm^/2. The dotted contributions 
involve higher derivative terms, e.g. (0"'"(— iV — 'V)/mB4'y which are suppressed by l/iVg. 
Again, note that the limit Nc and A large do not commute for a. For fixed A and large Nc 
the repulsion is strong as it should, while for fixed Ac and large A the repulsion is weak. 
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Figure 6: RPA contributions following from (165). 

Here V is the probing baryonic vector source. In large Nc the baryons are scalar fermions 
with no assigned spin to leading order in l/Nc- Therefore their spin degeneracy is 1. 
In the RPA approximation, the zero modes in (165) integrate to the effective action 

SRPAiV) = ^VlAlVl + ^VtAtVt , (166) 

with Jl = AlVl and Jt — AtVt for the longitudinal and transverse currents. ujAl^t are 
the longitudinal and transverse baryonic conductivities respectively. The RPA contributions 
as shown in Fig. 6 resum to 

A. = - 



1 - a^^UL ' 



At = , (167) 

1 - aKUr 



with 

nL,T = X]- +Ef^yAH^ + g)AHA;). (168) 
^-^ TUB ^-^ \ niB J 

k k ^ ^ 

The sohd hues in Fig. 6 lie in the Fermi surface, a/a 1/A^c- The transverse contributions 
follow solely from the dotted terms in (165). The summation is carried over the Fermi 
surface. is the massive and non-relativistic fermion propagator associated to (165) in 
the presence of a Fermi surface. The first contribution in (168) is from the seagull term in 
(165) and the second contribution is from the particle-hole bubble. The londitudinal vector 
response Al relates to the scalar density-density response function by current conservation. 
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Specifically, Hl = uJ^/q^H where 

F 

U{q)^J2^F{k + q)AF{k) (169) 

k 

is the Lindhard function for scalar fermions. For u!,q ^ but fixed /3 — uj/q/vp it takes the 
form 



n(g)«.^f-l + ^ln 



1 + 13 



1-/5 



The quasiparticle spectrum following from (169) is schematically displayed in Fig. 7 with 
massless quasiparticles of energy cu — vpq at small /c, and free massive fermions with uj — 
(f /2mB away from the Fermi surface. This description is rooted in weak coupling. 

We note that the longitudinal current in (167) develop a massless poles for strongly 
repulsive fermions as D4/D8. The longitudinal modes are directly tied to the Lindhard 
function by current conservation. They follow from 1 = all. In particular the longitudinal 
sound modes is stable for /3 > 1 above the quasiparticle cut (no imaginary part) with a speed 
iP » 1) 



cl ^ pvp ^ ]J ^'^^^^ vf^ V auB/mB = ci , (171) 

in agreement with the pressure p ^ an\/2 in leading order in the density. For (3 > 1 the 
decay of the longitudinal baryonic in D4/D8 is not through Landau-likc damping. 

The transverse mode follows from the pole at u;^ = aq^Hx- For a ~ a*/ {N^q^) the 
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transverse mode is gapped since 



Ub , riB 



ie{i-m 



SnuB 



(172) 



niB 20mB/3^ 



8ms ' 



for a;, g — > and fixed /3 = uj/q/vp- The transverse gap is typically ujt ~ a^nB/ruBNc with 
a* following from a a numerical analysis of the transverse quasi-normal modes in holography. 
In light of the RPA analysis, the holographic result suggests that the transverse baryonic 
fluctuations are long ranged and unscreened unlike their longitudinal counterparts. 

8 Conclusions 

We have analyzed the baryonic transport in D4/D8 (chiral) and D3/D7 (nonchiral) at finite 
density and/or temperature. D4/D8 is a holographic model of QCD at large Nc and large 
'tHooft coupling A. The transverse baryonic current in D4/D8 is saturated by the medium 
modified vector mesons in the confined phase with T < Mkk/t^- The vector spectrum 
is gapped since matter is uncomprcssiblc at large N^.- Confined D4/D8 matter becomes 
compressible to order 1/Nc with the occurence of a gapless longitudinal vector mode. While 
1/Nc effects are difficult to assess in holography, we have provided an RPA argument for the 
speed of the gapless mode using an effective action for baryons constrained by holography. 
D4/D8 is diffusive in the deconfined regime. 

D3/D7 is diffusive at all temperatures except zero where it is visco-elastic. This is a 
hallmark of holography at large Nc and large A. Indeed, the temperature effects are mediated 
by the BH background and leading in while the baryonic density effects are carried by 

the probe branes which are Nf/Nc suppressed. At strong coupling and large N^. the thermal 
or coUisional collision regime is dominant. The exception is D3/D7 at zero temperature but 
finite density as recently pointed by [17] . A Fermi surface (albeit strongly coupled) maybe at 
work in this case that suggests a visco-elastic regime. A longitudinal gapless mode emerges 
with a small width suggestive of a shear viscosity to baryon ratio tj/ub — h/A in cold but 
dense D3/D7. This mode is turned diffusive by arbitrarily small temperatures at strong 
coupling. Our observations extend readily to massive quarks in D3/D7. 
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A Cold Dp/Dq 



It is interesting to analyze the equation of state of cold Dp/Dq embeddings. We consider 
Dq probe branes whose worldvolume spans an AdSp+2 factor and wraps the n-sphere 5"" in 
AdS^ X S^, where p — q — n — 1. For example, q — 5,p — 2 corresponds to D5 branes on 
AdS4 X S^, and 5 = 3,p = 1 corresponds to D3 branes on AdS^ x 5"^ in AdS5 x S^. 

At zero temperature and for massless quarks, the pressure P and energy density e read 

[17] 

fl la ^2+1 

P = -TT = 7^^V , (173) 

_ _ pa _£±i 

e=-P + //BnB = -^^nB^ , (175) 



/ 1 X l/p / ^ £±1 



eo 47r3/2 \^r(p/2 + l)7riVeiV/y I 3^, 



(176) 



"9 



^o = 2y^('?^V''^nr^/^ (177) 



4Ar^Ar^ J p+1 



where = ^Ug, JIb = ^l^q^ and a = r(i/2-i/2p)r(i+i/2p) ^ ^ NjToqVn with K the 
volume of a unit n-sphere and T^g is the brane tension [18]. That is Nj = XNjNd {2nY, 
— ^^^3^ ; and A^3 = ^^^^ . The ratio e/eo shows the energy density in cold Dp/Dq 
normalized to the free Fermi energy density: 

- ^ X^'P . (178) 
Co 

For Dp/Dq = {D3/D7), {D2/D5), {D1/D3) dense Dp/Dq is unbound at large A. 



B Fermionic Drag 
B.l Hot D3/D7 

Hot D3/D7 is diffusive at all temperatures. The drag coefficient rju is inversely proportional 
to the baryonic diffusion constant Dg through the Einstein formulae at strong coupling [22, 
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23] 



VD = ^, (179) 

and for slowly moving particles. In equilibrium, the diffusion constant ties with the baryonic 
conductivity aq (at zero frequency and momentum) 

^_ <(A7V)^> _ a, 
TVs 

where S is the baryonic susceptibility. Thus 

Vd^-T . (181) 

The conductivity aq has been obtained using Ohm's law in [24] (See also next Appendix) 
^^^Mfl^^.^ d.^^, (182) 

with c = cos^^^(2;^,) and z^, = 4/(7r^T^) a dynamically generated value of a scalar profile at 
the BH horizon. Massless quarks correspond to 6* = and infinite mass quarks to 6* = n/2. 
Thus 

This expression expresses the drag of a quark in diffusive D3/D7 for arbitrary mass, tem- 
perature and baryon density. When niq = 

E 47r 27rT^ , 

where S = ^s^^^ [25]. When niq — oo 

^D = ^^ = -fVArS (185) 
where S can be read off from Eq.(5.7) in [25]. 



38 



B.2 Cold D3/D7 

The zero temperature case is visco-elastic as we suggested earlier. In this regime the fermionic 
conductivity ties to the diffusion constant by the Kubo formulae 

CTg = npDq , (186) 

+ (187) 

where / is Boltzmann distribution function and Ek — -y^W+m^. This relation follows from 
the relaxation time approximation in the quark probe phase space irrespective of strong or 
weak coupling. Relaxation to equilibrium at strong coupling is subsumed. For infinitesimal 
temperatures and for finite quark mass [24, 26] (see also next Appendix) 



ag _ NcNf 



T 47r 
Prom (179) it follows that the drag is 

rjD np 4:71 



Vc6 + d2 . (188) 



(189) 



This is the general form of the quark drag in a cold holographic medium (Coulomb phase) 
with infinitesimal temperature. We will assume it also for T = by continuity. 
For ruq finite (c 7^ 0) and Nc, A — > 00 (d — > 0) 



np 47r 



When T = and mg = 0(c = 1) 



where we used 



2(2-p) (l-p/2) 

VD = ^ , , (191) 

r(p/2) p-i 
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When T^O andruq^ 0(c 1) and p = 3 



rjoimg = 0) 



m.„ 



1 + 



a/ 



In 



(193) 



where = y + kp . 
For TTig — > oo(c = 0) 



47r riiT 7rVA^3 ^ 1 t\^T^ 



A^cA^/ d rig 2 2(p - 2) ' 



(194) 



where Ug = N^NfUB and 



rin 



For p — 3, rjD is the drag coefficient reported in [23]. 



(195) 



C Baryonic Conductivity 

The baryonic conductivity ag in D3/D7 has been derived by various methods [27, 24, 26]. 
Generically, the Kubo formulae for the conductivity is 

<7, = -hmilmG-*(i^) (196) 

a)->0 UJ u>=\k\ 

where only the transverse response function contributes, as the longitudinal part vanishes 
for light-like momenta by charge conservation. For a rotationally symmetric medium, G^^ — 
Gyy — Gzz are the components of the jxjxi jyjy a-^d jzjz retarded baryonic current correla- 
tions. 

Using AdS/CFT the transverse response can be extracted from (117), i.e. 

" »(» + d^»(-3 + 7»')) ; , tv'^u 1 + ud'' _ 
~ 2{1 - u^){l + d^u^) + (1 - 1 + u^d^ - ° ' ^^^^^ 

with u = . The horizon(-u = 1) is a regular singular point and the solution behaves as 
ot ~ (1 — m)^*™''^. We choose the incoming boundary condition(ar ~ (1 — m)^'"^/^) and 
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extract the singularity at w = 1 by substituting 



aT^{l-uY''°/^F{u). (198) 

F{u) is regular at u — 1 and satisfies the following equation 

„ u{-A + ud\3 - 7u^)) + 2i{l + m)(1 + dV)ro 

^ 2(l-M2)(l+d2M3) 

i(l + ^)(2 + dV(3 + 5^))ro + (-1 + ^ + dV(4 + 3^ + u^))to^ 

^ 4(l-^/2)(l+^,)(l + d2^/3) ^^^^^ 

In the hydrodynamic region(tti <^ 1) we may expand F{u) in terms of tti as 

F = Fo + mFi + m^Fs • • • , (200) 

and unwind Fq, Fi, F2, ■ ■ ■ order by order. For that consider the case with d = 0. The zeroth 
order equation is solved with 

^» = 1^ • P") 

where C is an integration constant. Regularity a,t u — 1 forces C — 0. So Fq is a constant. 
At next order we have 

^1' = ^^T-rV^^o + ■ (202) 

2(— l + K^j + 
Again regularity at = 1 sets the constant C = —iFo/2. Thus 

which is enough unwinding of the transverse solution for the Green's function in the zero 
frequency limit needed for the Kubo formulae. 

To obtain the retarded Green's function we need the boundary action 

lim -2N [ ^Ahh^a'j.{Z,K>)aT{Z,-n,) 

Z^oo J \2T[)^ XO 

- lim-4iV(7rr)^|^i^a^(K,ra)a^(«,-tt,) 
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where we recovered 27rQ;' and N = N{2T:a'Y = ^^ji- The retarded Green function is then 



G 



ret 



5^S 



5a^(0, tt))5a^(0, -m) (5F(0, ra)(5F(0, -m) 
^ ^ F(0,m) ' 



(204) 



and the conductivity is 



- hm -ImG'^";(X) 

a;— >0 UJ 



ui=\k\ 



87V(vrrrhmllmf^5(^ + 0(^- 
An 



(205) 



Extending the procedure to finite density yields 



NcNfT r- 



(206) 



where d = ^ = nJ^Vxt^ ^ ^^^^^ (^^2) ^^^^ ^ = ^• 

D D4/D8 Deconfined phase 

In this section we enforce an eigenmode analysis on the longitudinal and transverse vector 
currents in the deconfined hot and dense D4/D8. The eigenmode analysis parallels the 
one for cold and dense D4/D8 with reflective boundary conditions at the BH horizon. No 
imaginary parts arise from this analysis. In a way, in D4/D8 we may still entertain the 
possibility of stationary solutions between the Left-pending and Right-pending branes to 
mock up existing light bound states. Of course, this is a formal suggestion. 
The longitudinal operator (23^) is 



_ ^g^£T)A3 
= Oz—rTr-r, wzr^Oz + 



i^l/6 



IC^ 2 
K y 



(207) 
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When g = or a; = it is easily diagonalized, since 

2)L(g = 0) = -^dzy/KiK - l)Adz + , 

1 i^3/2 ^1/6 

= 0) = --a.^A»a. + ^A. (208) 

The Green's function (D^^) may be expanded in terms of the complete set of eigenvalues 
that diagonalize 

■^L{q = 0)/ = —y===A Xf , (209) 



= 0)/ = (^^^^A ) A/ , (210) 



where ^a^^rzi A and ^2^^^ A are weight factors. Using the complete sets, 



we have 



(Z|DZ^(g = 0)\Z') = Y: ""l^flf^l + ^^^^^^ , (211) 



(Z|DZ'(u, = 0)|g) ^ ^ ^"f + . (212) 

The transversal operator (Dt) is 

-Dt = idzVK(K-l)Adz + ^j==^ Uuj' - A-i-^gM \ (213) 
When 5 = it is easily diagonalized as 

■Driq = 0) = -^dzVK{K-l)Adz - £^ A . (214) 

V -f^ ~ 1 
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With the eigenf unctions and eigenvalues: 



{dzy/K{K-l)Adz)C 



,n 



the Green's function is expanded as 



{Z\D^\q^O)\Z')^Yl 



(215) 
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